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ABSTRACT 


This thesis offers a new and fundamental approach to the study of 
ship motion in a random seaway. The central idea is that the randomness 
of ship motion is caused by the forces and moments of wave action in the 
six degrees of freedom of the ship. These six forces constitute a multi- 
dimensional random process, which can be analyzed with recently developed 
mathematical theory. The advantages of this approach include: (1) The entire 
statistics of the sea as they affect ship motion can be determined from 
measurements made aboard ship, and (2) Theory is available for using the 
resuits of the analysis of these forces to develop optimum continuous motion 
predictions (as would be useful, for example, in an automatic landing system 
for aircraft carriers) as well as optimum control systems for the simulta- 
neous control of many stabilization power actuators, such as fins or rudders. 


In the thesis, mathematical models of the ship are presented in the 
Laplace transform and in the matrix differential equation forms. The pro- 
blems associated with correlation function and power spectral density appro- 
ximations are discussed and contemporary methods analyzed. An original 
design procedure is developed which will allow simultaneous maximum -effort 
control of n stabilization devices, minimizing a generalized quadratic error 
criterion, and which will take account of (1) all statistical interrelations 
between wave forces, and (2) the inter-coupled nature of ship motion. 


Finally, a comprehensive research program is proposed in the areas 
of ship behavior at sea and ship motion control which emphasizes the use 
of multi-dimensional random process theory as a basic tool, 
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Title: Assocıate Professor of Naval Engineering 
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CHAPTER I. 


ANEW APPROACH TO SHIP MOTION 


1.1 Introduction 


This thesis is concerned with a new approach to the understanding 
of the motions of ships in a seaway. A unified theory will be established 
which is heavily dependent on the methods used by the control systems 
engineer, who has spearheaded much of the contemporary advances in 


the study of complex dynamical systems. 


Complexity is the necessary penalty for an adequate analysis of 
the interdependent motions of a ship in six degrees of freedom. With the 
advent of high-speed, large capacity digital computing equipment, com- 
plexity and dimensionality alone are not formidable deterrents to under- 
standing, providing that the analyst has convenient methods for expressing 


the relations involved and has a well-defined objective for analysis. 


This work is to be considered as a practical application of general 
theory developed in a companion SURE. the latter being motivated by 
the dearth of theoretical aids to understanding of the random nature of total 
ship motion. Thus a major goal of this thesis is to interpret and to illus- 
trate the practical application of the mathematical concepts presented in 


Ref. 1, 


The hoped-for result of this work is to outline a new approach to 
ship motion which will yield dividends in more complete knowledge of the 
properties of the ship and of the sea, and what is perhaps just as important, 
will provide rational design procedures for control systems which tend to 


constrain or control ship motion. 


1.2 Current state-of-the-art 


A very useful summary of contemporary knowledge on the subject 





су 
of ship motion is found in the monumental work of Korvin-Kroukovsky , 


who has made a critical evaluation of most of the facets of past research 
on the sea and its effect on ships. It would be pointless to attempt to parrot 
this work to any extent, but there are many results of recent research 


which tend to support the approach to be presented in this volume. 
(a) The hypothesis of linearity 


First of all, consider the ship at sea, subject to hydrostatic and 
hydrodynamic forces which cause complex motions of the rigid structure. 
Conventionally, this motion is resolved into translations along three normal 
axes of a three-coordinate system and rotations about these axes. As will 
be illustrated in Section 2. 2, these motions can be considered as tbe solu- 
tion of a six-fold set of intercoupled second-order differential equations, 
with the sea providing the exciting forces. It is well-known that most 
systems which occur in nature can be approximated in a sufficiently small 
region of validity by linear differential equations, regardless of strong 
non-linearities in other regions. The question which naturally arises is, 
can the ship be adequately described by linear differential equations in its 


normal operating range of motion at sea? 


Des pite the prolific investigations of various non-linearities of 
ship motion which have appeared in the literature, it appears that the 
linear approximation is valid in many cases for the actual ship at sea. To 
support this hypothesis, model investigations made by Korvin-Kroukovsky 
and | and later by Be in regular waves showed good agree- 
ment between experiment and linear theory. cui Lewis and Nurn cad 
and Lewis and Ра1ге11' conducted towing tank tests in irregular seas which 
verified the principle of superposition, the consequence of linearity, for 
the motion of the model with respect to wave components at various fre- 
quencies. Also, in Section 3.2 of this work, statistical properties of mo- 
tion measured aboard the USS GYATT (DDG 1) in February 1959 by Marks 


and Dirkovie will be shown to agree nicely with a linear hypothesis. 





As an aside, it is indeed fortunate that the blessings of linearity 
can be attributed to the ship at sea, for a rigorous theory for non-!inear 
analysis of complex dynamic devices is in its infancy, and when stochastic 
processes are allowed as an excitation force, the understanding of such a 


problem is many years away. 
(b) Description of the sea and ship motion as a random process 


The sea defies a determinate description. Its ever-changing surface 
with waves rushing here and there, combining, disappearing -- this would 
appear to be a phenomenon which is absolutely indescribable in the large, 
However, some progress has been made in delineating the nature of the 


sea with statistical methods, 


The fundamental building block for contemporary statistical study 
of the sea is the unidirectional regular wave. If Superposition is assumed 
to apply, which regards the air-water interface as a linear transmission 
medium, the sea can be considered to be constructed of a distribution or 
spectrum of waves of various amplitudes, wavelengths, phases, and direc- 
tions, Since the sine wave has the virtue of orthogonality, has been conven- 
tionally used in the description of random processes in other fields, and 
is an approximation to the waves which propagate naturally at sea, it na- 


turaliy is used as the fundamental wave in the superposition procedure. 


Some of the major contributors to this ever-expanding body of know- 
ledge have been St. Denis and Person Neumann Denby bun MN Voz- 


1 13 
nessensky and Firsoff 2 Barber , 2. апа Marker 


е анаа ие has made what is perhaps the most appealing 
derivation of the description of the sea. He shows that a generalized auto- 


correlation function of the sea can be expressed by 


t 
(x, y, t) = hm 1 № (х, у, t) Ne y Ty, t +t) 


\ 5. ХҮТ | 
X, Y, To mom dx' dy' dt' (1.1) 





Here, N is the wave height at a point, x and y are normal position 


coordinates, and t is a time variable. 


Unfortunately, there is no known theory which would permit mani- 
pulation of the above expression or its possible transform mates (spectra) 
-^ provided that they could even be measured -- in order to reproduce use- 


ful results in the study and control of ship motion. 


Ne has had a considerable influence in the theoretical and 
practical development of sea spectral theory. By combining a shrewd 
analysis of early experimental data with the dynamics of waves at sea, he 
developed an expression for the spectrum (or frequency content) of the 


wave-height at a point 
ES 
к, ЖЕТІ; 
бо) - - е 7 722 
رر‎ 
where the К. terms indicate constants. 
ot. Denis and Deren hypothesized a directional sea spectra of 


K _ Ко 
М (д, @) = — е --2 cos e 
-w6 W 
where © represents the direction of wave travel as measured from a 


reference direction of maximum wave travel, 


They proposed that if measurements were made on a ship model 
to determine the response of the ship, both amplitude and phase, toa 
si:nusoidal wave excitation which varied over a range of frequencies and 
a range of oblique directions to the direction of motion, the response of 
the ship could be calculated by weighting each response function by the 
amount of excitation at that frequency and from that direction, and summ- 
ing the total. This idea has had a mjaor effect in shapıng current ideas 


about ship response in irregular waves. 


An alternate approach to the problem of defining wave spectra at 
{ ; 1 
a point was presented by Voznessensky and Firsoff e They made ex- 


perimental measurements of wave auto-correlation functions, which are 


cm 





the inverse Laplace transforms of wave spectra, and found that they could 
all be approximated by an exponentially-damped cosine wave, whıch will 


be later shown to agree nicely with a linear hypothesis. 


Finally, it is necessary to consider how the random process of 
wave motion arose in the first place. A number of investigations have been 
made into the mechanism of energy transfer from wind gusts to wave motion, 
According to Korvin=Kroukovsky*, this problem is essentially unsolved, 
From basic considerations, however, the wave action that a ship sees is 
generally a result of random wind action over great areas of the earth's 


surface and over long periods of time. 
(c) Summary 


It is useful to summarize at this point the major facts from contemp- 


orary theory which will be of use in the presentation to follow, 


(1) The ship is in general suitably described by the assumption of 


linearity. 


(2) The equations most used for "theoretical" wave spectra are 


actually based on an analysis of early experimental data. 


(3) The usual approach in expressing motion of a ship in a random 
seaway is to consider it as the effect of a wave spectra exciting a linear 
system whose continuum of transfer functions can be determined by ex- 


periment from ship model response to oblique waves. 


(4) Correlation functions of wave records show a tendency to appro- 


ximate an exponentially decaying cosine wave. 


(5) The actual cause of the waves observed by a ship are a very 
large number of essentially independent wind gusts, the effect of each 
making a very small contribution to the wave height observed at a single 


point, 





1.3 Outline of a new approach to the ship motion problem 


The heart of the contemporary approach to ship motion in irregular 
seas is presented in Figure 1,1, Here the wind is shown as a disturbing 
variable acting on the sea to create a wave spectrum. The wave spectrum, 
in turn, acts on a ship to produce motion. The original cause, the wind, 


is such a random phenomenon and is distributed over such a large area 


WAVE 
WIND MOTION 
SEA SPECTRUM “мү. ре ME 
| 


Fig. 1.1 Usual approach to ship motion analysis 


that any attempt to measure its effect would be ludicrous, Therefore, in 
the large, the system or mechanism which converts the total wind disturb- 
ance to a sea spectrum is forever unknowable. The spectrum itself is 
capable of measurement, but the oniy full scale attempt to measure it at 
sea, Project swop!. was a very drawn out, complicated, and expensive 
operation to obtain the statistical properties of a small patch of ocean at a 
certain instant of time. The continuum of transfer functions which relate 
the sinusoidal response of the ship to an oblique wave can be measured 

in the towing tank, but this would involve a very elaborate testing program 
because of two variables involved, direction and frequency. It appears 
that the study of ship motion, by focussing attention on the causal action 

of a sea spectrum which is in practice unmeasurabie, has painted itself 
into a corner. This corner is the stagnation which must necessarily result 


in a theory which depends on parameters which are so difficult to establish, 


suppose that a return to the most fundamental principals is made. 
The motion of a ship at sea is mathematically approximated by a system 
of constant-coefficient differential equations. In these equations appear a 
forcing function of six dimensions -- the forces along and moments about 
the three coordinate axes. These forces represent the effects of the random 


seaway to the ship. Fig. 1.2 shows the sequence of actions considered. 





WAVE 
SZ SPECTRU SHIP FORCES MOTION 


Fig. 1,2 A new approach to ship motion analysis 


Let us assume, for the moment, that the equation coefficients of 
a particular ship are known. This had been the goal of a great deal of 
theoretical and experimental work in the past decade, and at this writing, 
it appears that a reasonable approximation could be made on the basis of 
theoretical analysis and model experiments for a given ship. Knowing these 
coefficients, it would be entirely within the scope of current technology to 
measure the displacement, velocity, and acceleration variables of the 
various degrees of freedom of the ship, weigh these numerically by their 
respective coefficients in the differential equations, and compute the effec- 
tive forces on the ship. That is, we would continuously produce on a ship 
at sea a record of the approximate six force variables which would be 
composed of suitable weighted values of measured ship motion. This six- 
fold record would be a random process, a so-called multi - dimensional 
random process. Providing that our approximations to the equations of 
the ship are valid, we would have a record of the total random variation 
of the sea as it has affected the ship over the interval of time considered . 
It is important to note that the ship "system" considered here includes 
the hydrostatic and hydrodynamic forces which would act on the ship if it 
were suddenly released from an out-of-equilibrium condition in an other- 
wise smooth sea. The random force inputs, then, are the integrated super- 


imposed effects of the wave system acting on the ship. 


It is contended therefore, and this is the most important statement 
in this work, that the fundamental goal of the analysis of ship motion at 
sea should be to investigate the statistical nature of these wave forces on 


the ship. 


The logical question for the reader to ask at this point is -- How does 


one analyze a record of six simultaneous inter-related random processes 





to extract useful information from it? A theoretical answer to this question 
1 
was presented by the author in a companion document which was essen- 


tially motivated by the particular problem of analyzing ship motion at sea. 


Some of the advantages to be expected from pursuing this line of 


investigation are: 
(1) A catalogue of sea descriptions can be built up for a ship class. 
(2) All statistical data used can be collected aboard ship. 


(3) A quantitative means for determining the effect of hull coefficient 


changes on ship motion in typical seas will result. 


(4) A means of continuously making an optimum short-time predic- 


tion of any or all components of future motion is provided. 


(5) Based on this approach, a design procedure is possible which 
yields near-optimum control of steering and stabilization devices, allowing 
(1) for saturation of the force transducers, (2) for the intercoupling of mo- 
tions in the ship, (3) for the statistical relations among the various wave- 
produced forces and moments, and (4) for the possibility of noise-corrupted 


information from the motion sensors, 


The following chapters will develop and elaborate on the theory and 
practical use of this new approach to the understanding and control of ship- 


board motion, 





CHAPTER II. 


MATHEMATICAL MODELS FOR THE SEA AND SHIP 


2.1 Mathematical tools required 


The study of such a complex phenomenon as ship motion at sea 
is not suitably undertaken with horse-and-buggy mathematical methods. In 
fact, the success of such an analysis will be highly dependent on methods 


which can convey large ideas with a minimum of notational difficulty. 


The first major problem which must be resolved is the most efficient 
handling of problems involving the dynamic behavior of linear systems. 
Experience in the fields of control and communications engineering has 
shown that the most useful expression for linear systems analysis is the 
transfer function in terms of the Laplace transform variable, s. Its use 
avoids the unpleasant task of convolution when the response of a system is 
to be determined from knowledge of its excitation. Suppose a transfer func- 
tion G(s) of a system is given by 


Z(s) 


G(s) = P(s) 





where Z(s) and P(s) are polynomials in s. G(s) is the ratio of the Laplace 
transform of the output to that of the input. The easiest way of visualizing 
the meaning of this expression is to consider the differential equation which 
effectively relates the input x and output y of the system. If ЕЕ is replaced 
by the operator 5", then this differential equation will have the form 


Р(5) у = Z(s) x 


Formally, the definition of the Laplace transform G(s) of a function 


of time g(t) is given by 
G(s) = Јо е 231 dt (2. 1) 


Throughout the literature of naval architecture and, in particular, 





ship motion, the "jw" notation is used in place of the fundamental transform 
variable, s, to designate the transfer function of a physical system. Super- 
ficially, this would seem to be a difference of small importance, but in 
practice it has demonstrated the extreme preoccupation of contemporary 
theorists with the ideas of simusoidal response of systems, a point of 


view which it is hoped will suffer in the treatment to follow. 


The second major mathematical tool which must be brought into 
play is matrix theory, often called the algebra of higher mathematics. The 
use of matrices permits handling, at least conceptually, dynamic linear 
problems of arbitrary complexity. With six degrees of freedom of the 


ship, matrix methods are a necessity. 


It is assumed in the following that the reader is familiar with the 


basic elements of Laplace transform and matrix theory. 


2.2 Equations of motion of the ship 


In this section the equations which govern the rigid-body motion of 
a ship in its six degrees of freedom will be reviewed, and manipulated 
into a form most suitable for analysis. The conventional derivations of 
the equations of motion will not be presented here, as they are available 


17 И 
elsewhere , but some motivation for their appearance will be given. 


By equating the vector force F on a ship to the time rate of change 


of momentum, 





MJ Ame) 
ъ= کے‎ a 
dt dt | V | 
and the vector moment Q to the time rate of change of angular momentum 


— d — 4 = — — 
я | 
a general set of equations results which gives the three-dimensional motion 


—— 
of a free body under the action of external forces. Here, V is the instantaneous 


velocity of the ship, Jl is the instantaneous angular velocity, and IM is the 


10% 





angular momentum vector. 


A suitable right-handed orthogonal coordinate system is selected: 
(1) Under equilibrium conditions at constant speed U x is fixed 


in the ship from the center of gravity in the direction of forward motion, 


yis to starboard, and z is down. 


(2) The angular displacements -- ® (roll), © (pitch), and W (yaw) -- 
are measured in a clockwise sense from the equilibrium positions about 


the x, y, and z axes, respectively. 


(3) Linear velocities along and angular velocities about the x, y, 


апа z axes are denoted by u, v, w and p, q, r, respectively. 


With this choice of coordinates, six equations result which relate 


the force or moment in a given direction to the resulting motion: 


FY = m (x + qw - vr) 
F = m(ytru- pw 
Е, = m(z+pv - qu) 
М. = I, P + qr (I, - 1) 
= 1, 4+ Fp (I, - 1) 
M = I, r + pq (I, -1) 


These equations are linearized about the equilibrium point, yielding 


К = mx 
= 7. 
F с пут ғ 
ol E O 
F = mz-mU q 
2 О 
М, = Гр 
M = тб 
E S 
M =I r 
Z 2 


where ПА is the mean forward speed of the ship. The other terms involving 


the product of two velocities disappear because of their zero equilibrium 


values. 


Next, the hydrodynamic and hydrostatic forces exerted by the mo- 


are 





tion of the ship are examined. Initially, on a body of arbitrary shape it 

can be assumed that accelerations, velocities, and displacements in or 
about any axis cause forces in or about any other axes. This would indicate 
that 108 coefficients in general would be needed to specify the system, a 
staggering thought. A considerable simplification results in considering 


the real ship: 


First of all, linear displacements in the x or y directions could 
not be expected to result in a force on the ship. Next, the symmetry of 
the usual ship form about its centerline plane prevents small motion in 
the vertical plane (heave, pitch, and surge) from inducing forces in other 
degrees of freedom. Also, the surge or motion in the direction of mean 
forward velocity can sensibly be regarded as being decoupled from other 


motions. 


The linear differential equations which result from these considera- 

tions are certainly approximations to the actual case -- the hydrostatic 

and hydrodynamic forces are in general non-linear, the coordinate system 
used is valid only in a certain range near equilibrium, and the assumption 
that the wave excitation forces occur independently of the motion of the 
ship is certainly debatable. With all these misgivings, we still plunge 
ahead and trust that the resulting theory will provide (1) useful results 

for a ship at sea and (2) a foundation for a more complete theory as in- 


adequacies are revealed by future experimental results. 


The equations of motion are conveniently expressed in Laplace 


transform notation, using matrices, 


H(s) vs) ғ. (8) : Fis (3), (2.2) 


Here H(s) is a 6x6 matrix which expresses the dynamic forces 
which act on the ship as a result of its motion. V(s) is a six-dimensional 
vector of the displacements of the ship from equilibrium. F (s)|and F(s) | 
indicate the forces of contrdled elements and waves, respectively, in each 


of the degrees of freedom. 


9 = 





To structure these equations, the following components make up 
the displacement vector v. They are grouped in the order given in order 


to emphasize the interacting motions. 


Vi --Roll, o У. 
Pitch, © ; and v 


2 --Yaw, Y ; Уз --Sway, y; N ee Heave, Z; c Be 


6 =- Surge, xX. 
H(s) is constructed from three sets of intercoupling motions -- 


H,(s), roll-yaw-sway; H,(s), heave-pitch; and H.(s), surge. 


Н.С с i O 
-“Ір---1---- 
оо 
O num | Ну) (2.3) 
s (Kj ЕТЕ ик Кф SK + skr $». + SKr 
H,(s) = S^Np " s Nọ SAN +T,)+ SNr Ny SN s Ny 
S Yp + s(Yp +m Uo) Ye +5, +Уф SY +m) + SY 
(2.4) 
me) = ы + 5 و + ر2‎ S“ Zà + So - ^ Us ) 
: SC eS Moe Hz s*(Ma + Ty) + SMa guum 
(2.5) 
2 
H.(s) = E (m + Хы.) T 3 kel (2. 6) 


The components of forces or moments in Е. апа F follow the 


order of V above. 


This set of equations is believed to be one possible approximate 
representation of the linear model for the ship, but additions or simpli- 
fications will not affect the general validity of the theory to follow in future 
sections. The coefficients in the elements are expressed in standard nota- 
от, with the important exception of a reversal in sign. It is strongly 
felt that general coefficients of motion used in differential equations 
should be positive, reserving the negative sign as an indicator of such 
particular phenomenon as negative damping or a negative spring constant, 


which are usual sign-posts of instability. 
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solving Eq. 2.2 for the motion of the ship, 


vis) = H s) (ж, (s)|+ F is) (2.7) 





where 1 | | 
Eo. [NUS | ше 
0 | H. (s)! 0 (2. 8) 
2 
Aucun ee 
| | -1 
0 | 0 | H, (s) 


Figure 2.1 shows the ship as a ''system'', being excited by wave 


and control forces. 


SHIP 
ES + y | V 
M 
WAVE Force T SHIP MOTION 
CONTROL 
e FoRCE 


Fig. 2.2 General model for multi-dimensional ship motion 


Besides the transform approach, there is an alternate way of 
characterizing linear systems which is appealing in its generality and 
compactness of notation. This is by means of the so-called matrix differ- 
ential equation, as found in, for example, Beer In this method, the 
primary quantity of interest is not the input or output, but rather the state 
variables of the system. These are a set of quantities which, if specified 
as initial conditions, would be sufficient to determine tne entire free re- 


sponse of the system as it settled to equilibrium. 


To motivate this approach, consider the simple first order differ- 


ential equation 


a 
dt 


where a and d are constants, and c(t) is an independent variable. If c(t) 


x(t) = ax(t) + d c(t) 


is zero, and x(0) is specified, the solution of the equation is 


x(t) = x(0) ей” (020) 


As 





Now, suppose we write the same equation in matrix notation. 


E A = ооо (2.9) 


where A and D are matrices with numerical elements, and x(t) and c(t) 
are state and control vectors, respectively. If c(t) is zero, and x(0) is 


known, the solution of this equation is 


та” a) (2. 10) 


| At. 
from the theory of matrix differential equations. e у is an nxn matrix, 


known as the matrix exponential, which is defined by 


=> 
З = 2 A+ E 
A=O 


i 





i ! 


The ability to express the initial condition response in such a simple 


form for linear systems of arbitrary complexity is a considerable ad- 


vantage to the analyst. 


The ship equations of motion will now be cast into the matrix differ- 
ential form, and the numerical elements of the matrices will be identified 
as functions of the ship coefficients. Also, the state variables of the ship 
will be defined. In the following sections the concept of state and initial 


condition decay will play a prominent part. 


From Eq. 2.2, 


Hs)  Vís)| - Е (s)| + Rus) 


(2.2) 
H(s) can be expressed as a sum of three matrices 
2 
Hs) = s M * s В + R (2.11) 
where 
Kotil K> К. 0 0 0 
pex V 
M - Ne Net 1, М. 0 0 0 Еа 
ү. Y. Y.+m 0 0 0 
р У 
0 0 0 Zjtm i 0 
0 0 0 Му М.%1 0 
а И 
0 0 0 0 0 X tm 
ü 





p r У 
A Ne 0 0 0 
В = Р 
Y +m 0% Y. x 0 0 0 
(213) 
0 0 0 2 Z -m U 0 
W q O 
0 0 0 M M 0 
W q 
0 0 0 0 0 Xol 
R= Kg 0 0 0 0 0 
0 Ny 0 0 0 0 
0 Ур 0 0 0 0 (2.14) 
0 0 0 2, Zag 0 
0 0 0 M, Mg 0 
0 0 0 0 0 0 


Rearranging Eq. 2.11 and substituting the expression for H(s) into 


Br 2.2, ^M V(s)] = - В v(s)| - R V(s)] HF (s)| чл (s) | 
с 


Multiplying each sideby M 
2 -1 -1 2T -1 
S V(s) = -sM B V(s) -M К V(s)]+ M Е (s) M Fy (s) (2. 15) 


The following state variables are defined: 


x, = Roll velocity, p Xo 7 ourge velocity, u 
x. - Yaw velocity, r x, - Roll displacement, $ 
2 | 7 

х = oway velocity, v хз = Yaw displacement, Y 

x, = Heave velocity, w Xo = Heave displacement, z 

х = Pitch velocity, q ху Pitch displacement, Ө 
The following identifications can be made immediately. 

s? ws) - s X, | і-1,2,...6 

5 V(s)| = X. | j| 1242) 085 46 


Considering the product RV(s)| it is seen that this is equivalent to 


E 





N 
e 
о 
о 
© 
< 


1 
0 Ny 0 0 V 
в ү) = u 
0 Yy 0 0 І "REE j 
0 0 2, Ze| V, j= 7, 8, 9, 10 
0 0 M, Ме 
0 0 0 0 


Хх. Ay | 
5 Хр - ә | 

Хо Ay 

d 25 


Assembling these equations into a single matrix equation, with 0 


indicating a matrix of all zero elements, 


\ 
а ]- MB | -MK Mm! м”! 
e 27, АН + fey. fe 
100000 | x bee Fy 
010000: т A ^ 
000100 1 = 
000010 ! 
(2.16) 
ог, in general form, 
UN x= Ax + DF + D F (2.17) 
dt ec ta N w 


with the numerical matrices A, D, and D,, being defined above as a 


function of the coefficients of the equations of motion. 


To summarize, the ship has been described in terms of a Laplace 
transform matrix, H(s), and a numerical matrix A-- the two basic des- 


criptions of linear systems --in terms of the physical parameters and 


E 





hydrodynamic coefficients of ship motion. This will permit the efficient 


discussion of ship motion in the following sections. 


2.3 Important relations in scalar random process theory 


The purpose of this section will be to present and analyze some 
pertinent relations in the study of random processes and their effects on 
linear systems. The emphasis here will be on viewpoint, and the reader 


1 
is referred to standard texts Y or Ref. 1 for proofs and extensions, 


Statistical analyses of random processes are made in order to de- 
tect relations which exist between signals at separated instants of time. 
The class of random processes considered here are stationary -- the 
statistics do not vary with time -- and satisfy the ergodic hypothesis, 
which essentially means that the desired statistics can be adequately 
approximated through analysis of signals over a long interval of time. 

This theory becomes more applicable as the process considered approaches 
a Gaussian random process -- which latter case is approximately true in 


the ship motion problem. 


The fundamental measure of random processes is the correlation 


function, (p (T), defined by 


и = Е EG . y(t m (2. 18) 


where x(t) and y(t) are two not necessarily distinct random processes, 

and E B indicates the expected or mean value. When x and y are different, 
Y en is known as a cross-correlation function, and when y is identical 

to x, ф (т) is an auto-correlation function. Fig. 2.2 shows typical correla- 


tion functions. 


Pe er) (>) ty т) 


T » 


Fig.2.2 (a) Typical auto-correlation function. (b) Typical cross-correlation 
function. 
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Since correlation functions (for T70) have the appearance of transient 
signals in a linear system -- that is, decaying to zero as F— oO --it 
is natural to use their Laplace transform as an alternate measure of согге!а- 
tion, These transforms are known as auto or cross power density spectra. 
Since correlation functions are defined for positive and negative time, the 


bilateral or two-sided Laplace transform is conventionally ee. 


If jo is substituted for the Laplace transform variable s, the auto 
power density spectrum can have the connotation of an amplitude-squared 
distribution of sinusoidal frequency content in a signal except for a constant. 
For a cross power density spectrum in jo, the information is conventionally 
represented in naval architecture literature by co and quadrature spectra, 


which are its real and imaginary parts, respectively. 


It is the contention of the present author that preoccupation with 
the frequency or w representation of random processes tends to stifle 
progress in the understanding of random processes, especially when cross 
power density spectra are considered. For example, it is impossible to 
acquire much feel for what is actually happening between a pair of random 
processes from consideration of their co and quadrature spectra, yet these 


observations have appeared quite frequently in recent literature. 


As amore basic statement of the worth of power density spectra 
in the analysis of random processes, this thesis will demonstrate that 
their main use should be as a stepping-stone in determining the transfer 


functions of linear systems. 


The general formula which governs the change in statistical prop- 


erties of signals after passage through a linear system is 


G(s) = 2 2GCs9H(G Y xy, (s) (2. 19) 


Az! j=1 
Here, Q. (5) is a cross power density spectrum between signals 
x and y (they may be identical) or the Laplace transform of PD: xis 


made up of n signals, х., each passing through a linear system with trans- 
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fer function G(s), and y is similarly composed of m У; signals, each of 


which have been operated on by H.(s). The set of input statistics 
Ф хуб) т) 
is assumed known in order to completely specify x and y. 


As a trivial example of the application of this formula, consider 
Fig. 2.3, where the input spectra of x, Ф. (8), is specified, and the out- 
put spectra Ф в) and the cross spectra between x and y, È (s), аге 


to be found. 


Н (5) 


Fig. 2.3 A simple linear system 


Since ху = х, G, =]: E and H,(s) = H(s), 


Ф (s) 


IY 


H(-s) H(s) Bs) 


and 


H 


(s) 


=) 


His) Ф (в) 


To make an application of this result which will be of continuing 
Significance in future sections, suppose that the input signal x has an 
auto-correlatiqn function which is a unit impulse at the origin. Since 
Ф (0) = E (oy, this impulse would mean an infinite mean-square 
value which is, of course, not possible. However for many practical 
problems the assumption of an impulse is valid, just as in mechanical 
problems a force impulse is used to approximate а delivery of momentum 
over a small finite interval. This signal is known as white noise, and in 
the transfer domain, has a constant value at all values of w . It is observed, 
for example, as ''grass'' on a radar scope. The important thing about white 
noise is that is the ultimate in randomness and is totally uncorrelated 
with time-shifted versions of itself. Also, it can be approximated very 


nicely for use in analog simulation. 
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Thus if © (s) = 1, then 
XX 
® (s) = H(-s) H(s) 
Dy 


It is known that power density spectra of single random processes 


expressed as rational functions of s can always be factored into a form 
È (s) = H(-s) H(s) 


Thus, a random process can always be viewed as the output of a 
linear system, denoted by H(s), which is excited by white noise. By postu- 
lating the totally random excitation, one can focus attention on the effective 
linear system rather than on the statistics of the random processes. This 
viewpoint is not just an abstraction, for an hypothesis has been offered 
in the ur that most random phenomena observed can be traced 


back to an original cause which has the properties of almost-white noise. 


An appreciation for the strength of this approach can be obtained 
through considering two fundamental problems of random process theory 
--the optimum filter and the optimum predictor. These new results were 
first presented in Ref. 1, and illuminate the basic mathematical solutions 


given in the classical work of ia 


a. The optimum filter problem 


Figure 2.4 shows a linear system W(s) which must operate ona 
random process v. The signal v is assumed to consist of signal s and 


noise n, with the criterion for optimum performance being the minimiza- 





Fig. 2.4 The basic filtering problem 


tion of the mean square error between output and signal, e* 
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If noise and signal are uncorrelated, 


Ф (з) - © (в) + © (5) 


from Eq. 2.19. Let us assume that the factors of the denominators of 


each term the poles, are distinct. Suppose È (s) is factored such that 
Ф (з) = G(-s) G(s) 
VV 


with the poles and zeros of G(s) all located in the left half of the complex 


plane. 
Also, let G(s) be separated into two terms 
G(s) = S(s) + N(s) (2.20) 


where S(s) contains all the signal poles and N(s) contains all the noise 
poles, that is, the positive poles of the signal and noise spectra. The ran- 
dom process v could be considered to be constructed from white noise w 


as shown in Fig. 2.5. 





Fig. 2.5 Effective generation of noisy input random process 


Ref. 1 proves that the optimum linear operation on the random 
process v merely reproduces the output S of Fig. 2.5 as the best estimate 
of the actual signal s. This optimum system is shown as a simple unity 
feedback configuration in Fig. 2.6. This representation is also valid when 


signal and noise are correlated. 





Fig. 2.6. Simple feedback configuration for optimum filtering 
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This is the first example of the utility of postulating a model for 
creation of the random process, contrasting with the laborious mathe- 


matical treatment found in standard texts. 


b. The optimum predictor problem 


It is desired to produce continuously the best estimate of a future 
value of a random process in a mean square sense. This problem (more 
specifically the prediction of aircraft motion for fire control applications) 
was the motivation for Wiener's wartime ver which laid the foundation 


for modern random process theory. 
If one postulates a model G(s) for the process v, where 
G(-s) G(s) = Ф, (5) 


then this transfer function model is characterized by a state just as any 
linear dynamic system is. These state variables are analogous to a suffi- 
cient set of initial conditions needed to determine the future output at v 


if the white noise excitation were suddenly cut off. 


The optimum predictor, as proved in Ref. 1, is one which con- 
tinuously produces or recovers the state variables of the model G(s) and 
weights each by its value after T seconds of initial condition decay. What 
this means is that the best guess of the future value of a random process 
is made by neglecting any future effect of white noise and only by con- 
sidering the future action of present signal levels in the assumed model 
G(s), which summarize the total effect of the past theory of the random 


process as it can affect the future. 


The importance of this result transcends the specific application 
of a predictor. Since the expected value of future signal in a random pro- 
cess is found by allowing the random process state variables to decay as 
initial conditions, these quantities assume prime significance in the design 
of a control system to operate in the random environment. These matters 


will be considered in some detail in Chapter 4. 
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2.4 Multi-dimensional random processes 


By far the majority of theoretical investigations into random pro- 
cesses and the optimum systems which operate on them have been con- 
cerned with single or scalar signals. When two or, in general, n random 
processes are to be processed simultaneously, the theoretical aids to 
understanding have been essentially missing. One of the major purposes 
of Ref. 1 was to meet this need, especially motivated by the desire to examine 
random behavior of ships at sea. This section will outline the problem and 


a general theory for multi-dimensional random processes, 


First, it is necessary to consider what a multi-dimensional system 
is. We will concern ourselves with linear systems which have inputs and 
outputs of n signals. Fig. 2.6 shows the structure to be considered. Since 
each output а; сап be caused by any or all of the inputs ШЕ an entirely gen- 


eral representation of W(s) is the definition 


4 .(5) 2 W..(s) v.(s) O A a 
1 БЕЗ 1) J 


or 


4()- wt) v] (2.21) 


LD © © 
г = 


ee 





Im 


Fig. 2.6 A general multi-dimensional system. 


Each element of W(s) is a separate transfer function. Fig. 2.7 


demonstrates the configuration of a 2x2 matrix system. 


a) (мө мө | 
Ч, 7 Wa (5) Wals) V 





Fig. 2.7 System defined by 2x2 matrix 


Од = 





Next, we consider a single set (or ''vector'') of signals, v, and 
inquire as to what information must be known in order to completely specify 
it as a random process. Obviously, from the previous section, the only 
measures we can use involve correlation and thus specification of the 
power density spectrum Qv v (5) fori, j= 1,2. . . n is sufficient to 
define the process statistically for our purposes. In matrix form, an nxn 
matrix with ijth element Qv v (s), Q (=), represents v, Similarly, the 


croSs-power spectral relation Qv. ш) is the ijth element of a matrix Ora! 


A general formula for the matrix of cross-power density spectra, 


D (s), between two vector random processes, x and y, is 
— AA 





2% jr. 
E К T 
Ф. (5) = = 2 С. 5) QD x x (s) H. (s) (2929) 
where yx 
x(s)| = A G(s) X (в) | 
апа Y(s)| И 2 Ne) Y. (s) 
J) Ceo c J 


The superscript [ ie indicates the matrix transpose. This simple 


formula is analogous to the scalar case, Eq. 2.19. 


To illustrate a simple application of this result, consider again 


the system of Fig. 2.6. Here, 





б (в) = М(-) Ф (в) wis) (2.23) 
апа КЕ 
P. (s) =- © (s wis (2.24) 





туд AVN 
Suppose we next inquire as to the output relations of a white noise 
driven system, since this concept supplied so much insight into the scalar 
problem. Here, the white noise is an n-dimensional vector of uncorrelated 
white noise with unit power density spectra. D,, (5) in this case is I, the 


identity matrix. The output vector v of a matrix system G(s) is, from 


Eg. 2.23, 


MO = G(-s) E (8) (2. 25) 
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Let us now ask the reciprocal question --if ф (в) is found by 





measurements on v, can a G(s) be determined which satisfies the above 
relation and could create the process from white noise? This is a problem 
of matrix factorization and, once solved, it serves to unlock the door 


to the understanding of multi-dimensional random processes. 


A general method is presented in Chapter 3 of Ref. 1 which always 
leads to a solution to this problem. We will not concern ourselves here 
with the details, save to say that the procedure used is valid without regard 
to order and consists of elemental operations such as polynomial factoring 
and solutions of linear equations. However, the significance and utility 


of the answer is of considerable importance, 


The solutions to the multi-dimensional filtering and prediction 
problems are identical to those cited in the last section for the scalar 
case --1f matrices replace the single-dimensional system expressions. 
Especially important is the fact that the optimum prediction of future values 
of the multi-dimensional random process is the decay of state variables 


determined in the model G(s). 


2.9 Significance of multi-dimensional theory to ship motion study 


A foundation --admittedly sketchy --has now been laid for the 
presentation of a new approach to statistical understanding of the sea and 
ship motion. The heart of this approach is the visualization of a random 
process, single or multi-dimensional, as the output of some linear system 
which is driven by white noise. The randomness of the random process 
is caused by the white noise, and the memory or useful part is caused 
by the linear system. 

First, let us consider the sea itself. It is manifestly impractical 
to monitor a large area of ocean for amplitude disturbance, Project swop!? 
notwithstanding. However, a finite set of measuring instruments such as 


wave poles is possible. The signals from these devices comprise a multi- 
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dimensional random process and can be analyzed by cross-power spectra 
techniques to produce the matrix of spectra, Q (s). The G(s) which 
paa" a 


satisfies the equation 


GCs) G'(s) = © (в) (2.25) 





is the effective system which carries in its signal levels all that is know- 


able for prediction. 


White noise can be visualized as a random sequence of uncorrelated 
impulses. In the case of the sea, the white noise is physically approximated 
by the gusts or random wind forces which continually supply energy to the 
sea. The sea is then a linear system which converts the effect of a single 
gust into a small displacement disturbance which propagates at the inter- 
face and interacts with the contributions from millions of other gusts. This 
very gross description is intended to make the point that white noise is not 
necessarily an abstraction, and is actually a major factor in the actual 


randomness observed at sea, 


If we were to use an arbitrarily large number of wave poles, it is 
possible that we could reduce the prediction error for future wave motion 
at a given point to less than any desired amount, since the inclusion of 
an additional signal source which is correlated with the signal in question 


will always reduce the prediction error. 


These ideas are interesting but of dubious practical value since it 
is a rare physical problem which is concerned with instantaneous values 


of wave height. 


A second use of the multi-dimensional theory occurs when ship 
motion is to be measured. If the six-dimensional motion of the ship, q, 
is measured, and the input forces are denoted by v, then from Eq. 2.23, 


Q (9"H C9 DG n^" S 


-1 
where H (s) is the ship matrix transfer function as derived in Section 2.2. 
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Suppose first that H(s) is known, from knowledge of the coefficients 


of the linear differential equations. Then, the input spectra Ф ve) 15 





given by 
Q (з) = H(-s Ф, (5) н (в) 











An approach which is less interested in the mechanism of ship motion 


would take the motion statistics, Ф (9). and factor this such that 
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This system G(s), composed of a mixture of force and motion systems, 
has in its signal ET state variables the necessary information for 
the optimum prediction of future motion of the ship. This would be of 
considerable practical import for motion control devices, and for perhaps 
more exotic uses such as guided missile launching installations or auto- 


matic landing systems on aircraft carriers. 


However valid this parti cular approach might be for a certain spectrum, 
it is felt that the relatively constant transfer function of the ship should be 
determined and fixed in any design procedure, and attention be paid to the 
real variable in the problem, the random process of forces on the ship 


due to wave action. 


The ship coefficients are perhaps best determined in the towing 
tank, but a useful check can be made for some of the major parameters of 
interest --such as natural frequencies and damping ratios --by analysis 
of the correlation functions obtained from reduction of motion data taken 
on the actual ship at sea. Chapter 3 will illustrate an example of this par- 
ticular use of these statistical measures in motion analysis of the USS 


GYATT (DDG 1). 


The primary emphasis, then, should be first on fixing the para- 
meters which govern the differential equations of motion of the particular 
ship, and then investigating the statistical random process of the wave 


forces on the ship. The latter is best done practically by simulating aboard 
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ship the equations of motion. This is best shown by means of an example: 
Suppose that the roll moment is to be found, From Section 2. 2, 


М. = (к. 19% + К, ф ақ; TES ф +к,ф +K, J 


where ф, Y , and y are roll, yaw, and sway displacements, respectively. 
These motions can be measured quite easily with accelerometers, rate gyros, 
etc., and weighted and summed as shown. The control term M, refers to 
the moments generated about the roll axis from, for example, fins or rudders. 
These forces are probably best introduced into the problem from strain 
gauges on the rudder or fin stock or other signal sources which give an 
output proportional to actual lift force obtained, not to the deflection angle 


of the appendage. 


The result of this measuring and weighting for all dimensions would 
be a set of six random processes recorded simultaneously. The cross- 
power density spectrum Pr v (s) would be P between signals, and 


the resulting matrix Ф Ж factored into G(-s s)G en For the purposes 








of the ship, the seais is represented over the interval of 
data taking as a six-dimensional matrix filter G(s) excited by six independent 


white noise sources, 


A library of these G(s) representations could be established. Typical 
patterns could be determined and useful approximations could be made in 
their configurations as various sea conditions were "tested" by the ship. 
Hopefully, the sea could be well-characterized by a few parameters of 


the matrix. 


A very important use of the G(s) representation is that it is in a 
form suitable for direct instrumentation on an analog computer, using 
commonly-available white noise sources to excite this linear system. Thus 
the effects of, for example, a strong correlation between pitch and heave 


excitation could be experimentally analyzed for the purposes of control. 
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It is a reasonable assumption that G(s) representation of random 
seas would not differ greatly for hull forms of ordinary proportions in 
the same seaway. Thus, if a typical G(s) were specified after some ex- 
perience with this statistical measure of the sea, the ‘effect of coefficient 
changes of a particular ship could be evaluated through considering the 


matrix 
® aals) - H` (s) G(-s) Gl(s) n (в) (2.23) 


The diagonal elements of © (5) are the auto spectra of roll, pitch, 
etc. The mean square value of these quantities is the zero value of the 
corresponding inverse Laplace transform, Thus an analytical relation is 
available for studying the magnitudes of motion of a ship in a standard 


sea as a function of the coefficients in its differential equations. 


With the increasing use of digital computers for complex problems 
in naval architecture, it seems possible that quite an ambitious scheme 
for hull design could be formulated, From a set of equations for hull forms 
with several variables, the coefficients of the differential equations could 
be calculated by strip theory. These coefficients could be used to evaluate 
mean-square motion of the ship to the standard seaway G(s), and the result 
used inside a loop in the computer to find the optimum set of hull para- 


meters. 


The most important immediate use of this multi-dimensional ran- 
dom process appraoch is in control. This interesting problem will be dealt 


with in more detail in Chapter 4. 
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CHAPTER II. 


PRACTICAL PROCESSING OF RANDOM PROCESS STATISTICS 


3.1 Introduction 


In this chapter we will be concerned with the problem of obtaining 
useful information from actual measurements made on random processes 
in the form of correlation functions or their transform mates, the power 
density spectra. The theoretical results of this thesis, and indeed all ba- 
sic random process theory, depend upon the ability to approximate corre- 
lation functions by exponentials or damped sinusoids, or alternately, to 
represent spectra as meromorphic functions of s, that is, with numerator 
and denominator which are polynomial functions of s. This does not nec- 
essarily pose a significant restriction on the accuracy of the theory, how- 
ever. Laning and Batre for example, have shown that correlation func- 
tions can be approximated by the exponential class of functions to within 


any desired accuracy level. 


In advocating this approach, the author hopes to supplant the well- 
known Neumann spectrum, presented in Chapter 1, with a less rigid model 
for the presentation of wave-height statistics. In fact, the term e En which 
Neumann used in approximation of early experimental spectra, is analyt- 
ically untractable for the purposes of random process theory. The alter- 
nate results of Voznessensky and КОО © however fit in directly with 


this approach, using an exponentially damped cosine wave to approximate 


many measured wave-height auto-correlation functions. 


Also, in several examples of correlation functions of ship motion 
to be presented in this chapter, the exponentially-damped sinusoid will be 
seen as a fundamental building block. These particular statistics are re- 
presentative of measurements made aboard the USS GYATT (DDG 1) by 


Marks and DEES in early 1959. 
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3.2 Correlation function analysis 


In this section, methods will be examined for approximating a 


correlation function with the following terms 
р(9- 2C e i wp rt $) 


allowing the special case with В; апа ф. both zero to represent the pure 
exponential. Additional but less frequent allowable terms are of the form 
Ae S АН + $). This emphasis on decaying sinusoidal functions 

is made because it corresponds to the usual observed behavior of statistical 
measures of ship motion and wave-height records. 


The standard form for the transfer function of an under-damped 


second order system is —— — Вы . 5 is defined as the damp- 


2 
S + 250)8 i e the natural or undamped 


ing ratio, and 2: 
-at 

frequency. A typical transient of this system will contain the term e Ñ Cow Bt. 

a and ß are related to the standard parameters, J and wo» as follows: 


2 2 2 a 
к E u cs 


These relations are best visualized with the ‚aid of relations shown 
on a plot of the system poles on the complex plane, as shown in Fig. 3.1. 


S-PLANE 





Fig. 3.1 s-plot of standard second-order pole 


These standard parameters have been adopted because of their 
utility in plotting frequency response in control systems applications, The 
damping ratio governs the rate of decay of amplitude of oscillation. Fig. 3.2 


shows the amplitude ratio of successive positive and negative peaks as a 
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function of $ for the useful ranges of damping in ship motion. 


Another characteristic of the damped sinusoid, the distortion of 
the normal sinusoidal waveform because of exponential attenuation, is 


presented in Fig. 3.3. 


When it is desired to approximate a correlation function which has 
been experimenially determined by one or more damped cosine waves, 
the procedure becomes quite difficult if two nearly equal frequencies are 
involved, because of a ''beat'' phenomenon. In this case, the spectral or 
transform methods to be described in the next section become very attrac- 
tive. In the remainder of this section, some simple aids will be given to 
speed determination of the characteristics of single significant sinusoidal 


components of correlation functions, 


The first example to be considered is an auto-correlation function 
Dr (T) of the roll angle of the USS GYATT (DDG 1), measured in state 
five and six Computed values at one-second intervals are shown in 


Fig, 3,4, for positive values of T. 


First, the approximate times of peaks are determined and averaged, yield- 
ing 4.46 seconds as the mean time between adjacent positive and negative 


peak values. Since .46ß = т, Шеп В = . 705. 


The interpolated absolute magnitudes of %- T) at these peak times 
are nlotted on semi-log paper. Using the peak values eliminates the effect 
of the sinusoidal term; using a log presentation allows an exponential to 
plot as a straight line. Figure 3.5 shows these points and a superimposed 
straight-line approximation. The intersection corresponding to the time 
p of this line and the = ordinate is at time 17.1 seconds. Thus 


а = т = ,0583, The standard parameters are then J = . 082 s and 


O (ТҮ 
n 


E т 
The approximation e ' 222% сов .705T is plotted in Fig. 3.6, and 


Shows quite a good agreement with the measured values. It is heartening 
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to note that this auto-correlation function, which in this case is essentially 
a measure of the dynamic behavior of the ship in roll to a totally random 
excitation, corresponds to a linear theory even at large amplitudes of 
motion at sea. The peak amplitude of roll observed at the time these 
measurements were taken was 28.5 degrees, which is certainly out of 


the range of a small-signal theory. 


Fig. 3.7 presents a similar approximation to a pitch auto-correla- 
tion function, PN. The observed characteristics are У = .0869 and 
=. 052. The approximation is not as good as was that of the roll func- 
tion, which reflects (1) the non-linear aspects of the pitch system at large 


amplitudes and (2) an indication of the wave excitation statistics. 


Figures 3.8 and 3.9 show the measured cross-correlation between 
roll and pitch. Here do = Ф (7): а well-known identity . Because 
of the difficulty in simply expressing these function in terms of damped 
cosine waves, their analysis is deferred until the next section. However, 
it is important to note that the cross-correlation measured between these 
roll and pitch signals does not arise principally because of coupling between 
the roll and pitch degrees of freedom in the ship system, but rather from 
the cross-correlation between the exciting moments in roll and pitch from 


the same wave spectrum. 


The significance of correlation function analysis in the time domain 
is that dominant frequencies and damping ratios can be readily determined, 


as well as an indication of the non-linearity of the source of the signal. 


3.3 Analysis of power density spectra 


This section will discuss methods of approximating power density 
spectra by rational transforms in an efficient fashion. Also, conventional 
techniques for digital computation of spectra from correlation functions 


will be critically analyzed. 


The heart of this technique is based on methods used by electrical 
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engineers to approximate the frequency response of physical systems with 
transfer functions in the Laplace variable s. Consider first the auto power 
density spectrum, Q (s), which can be represented with the following 


factors in numerator and denominator: 
2 2 
К, (+з+а.), (в +2} 0). 5+ 2). ) 
ea г 


Suppose the logarithm of Ф (s) is taken. In illustration, if 
б ТЕ о (s+a) 
хх 2 2 


2 2 
(в -2% sto ) (5 t2Yo sto) 


log D , (s) - JO Me IE -s+ta)+tlog(s+ta) 
2 


2 2 2 
с0а (Бе 25 w St w A EA oS e ) 
In general, each separate pole or zero has an additive effect on 
the spectra when the latter is expressed in logarithmic form, and thus 


each can be considered separately. 


If s is set equal to jù, the mirror symmetry of poles and zeros 
about the imaginary axis in the s-plane for auto power density spectra 
ensures that each pair contributes a real value only. Thus the magnitude 
of each pole and zero as a function of 4 is the item of interest, and the 


phase angle is uniformly zero in this application. 


For quick graphical approximation, it is customary to plot the 
frequency function on a logarithmic frequency scale. Fig. 3.10 shows the 


ff f | | TS 
effect of a single simple pole, aa 
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Fig. 3.10 Frequency response of a single pole 
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Here, the asymptotic approximations reflect the fact that for small W , 


1 
GER А z 
Е . Because of the logar 
ithmic frequency scale, the high frequency asymptote plots as a straight 


Els = , and for large W, 
Jwra a 





line, which makes the entire frequency response of each basic building 
block shiftable in frequency without distorting its shape. As another example, 


Fig. 3.11 illustrates the magnitude approximation of a frequency function 
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Fig. 3.11. Approximation of frequency function with poles and zeros 


with one zero and two poles. Here the asymptotes are plotted all at once, 
and the positive or negative corrections are added to approximate the ob- 
served frequency function. It should be clear that the use of what is eff- 
ectively a log-log presentation has resulted in considerable simplification 


in a trial-and-error approximation procedure. 


The conventional unit for expressing the magnitude of the frequency 
function, F(jW) , is the decibel, defined as 20 log, 9 |F Go) . Although 
this unit has a meaning of power in audio applications, this restriction is 
not valid for this frequency response use and the given definition applies 
to any physical quantity described as a function of frequency. For an asymp- 
tote oí decay of a single pole, > , the slope must obviously be 20 db fall 


in a decade of frequency in W. 


Considering the case of auto power density spectra, the log mag- 
nitude of each LHP simple pole will be doubled to account for the pair of 
mirror-symmetric poles. A plot of magnitude of a simple pair of poles 


of a spectra is given in Fig. 3.12. The reminder is offered that multiplica- 
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tion of this spectra by a constant K corresponds to a upward translation 
of 20 108,0 K db, and that replacing s by = corresponds to a lateral shift 
of the plot so that the "breakpoint" --or intersection of the asymptotes -- 


occurs at "a'' radians/sec. 


Of more immediate interest for the analysis of wave height and 
ship motion spectra are the poles which are complex-conjugate and indicate 
sinusoidal behavior of the auto-correlation function. Fig. 3.13 presents 
the normalized curves for a spectrum of the form 

б К 1 


ГЕ я 2 =. «il eg 25622 +1 | 
n n n n 


over the important range of damping ratios for these applications, Note 
that the asyptotic slope is now 80 db/decade, indicating the presence of 


four poles. 


If zeros are to be used instead of poles, the same curves can be 


used with a reversal in sign. 


Before using these tools to examine Several sample spectra, it is 
necessary to review means used to determine these particular spectra 
from their associated correlation functions. The steps to be outlined are 
given in the authoritative work of Mare and are currently used at the 
David Taylor Model Basin to produce spectra from digitally computed 


correlation functions. 


A finite duration of auto-correlation function 0. T) is determined 


at the first M sample points, T seconds apart. Its Laplace transform 
ed 


ec) 
= | Be EE > 
far Y un Ze = COSWT -]j sin JT) = 2 far PT) COS WT 
is approximated by : 
hist 
aT zo (KT) cos W. kT (i=1,2...m) 
Et XX i 


This approximation is further smoothed by the operation 
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where the superscripts s and u indicate smoothed and unsmoothed respec- 


tively. 


An auxiliary purpose of the analysis of these spectra will then be 
to investigate the effectsof (1) approximating the infinite correlation func- 
tion with that of a finite length, (2) approximating a continuous integral 


by the trapezoidal rule, and (3) smoothing the rough spectra. 


The roll auto-correlation function of the USS GYATT (DDG 1), given 
in Figure 3.6, was excellently approximated by a damped cosine wave, 
which has a known spectrum. Fig. 3.14 shows the rough and smoothed 
spectral data computed with the above scheme, and the spectrum of the 
second-order poles with w = and | computed in Section 3.2. An excellent 
fit is obtained near the resonant frequency with the rough computed points 
of the spectrum. It should be noted at this critical resonant point the 
smoothed points indicate a greater damping than is actually the case. This 
15 the real penalty for such a crude smoothing operation and easily ex- 


plained. 


Fig. 3.15indicates two examples of this form of data smoothing. 
In (a), the smoothed points tend to approach the desired straight line. In 


(b), a peak in perfect data is considerably flattened because of smoothing. 


ROUGH SMOOTH 
X 
x 
a) 
v | 
O 
b) {\ 


Fig. 3.15 Effects of data smoothing 


With the resonant peaks expected in wave height and ship motion spectra, 
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the smoothing operation is useful in producing a continuous (if not true) 
curve away from the resonant point but causes a considerable loss of in- 


formation about the damping of the system under investigation. 


The points in Fig. 3.14 labelled ''N" indicate negative power density 
in the rough spectral computation, which of course, could not occur with- 
out imaginary signals. These negative points at low frequencies indicate 
the effect of a finite (30 second) length of auto-correlation function considered. 
At high frequencies, the approximation of the integral by a finite summation 
causes the rough spectral points to be alternately too positive and then ne- 


gative. 


To summarize, the representation of spectra digitally computed 
by contemporary methods can be considered essentially accurate only in 
the determination of resonant frequency, less accurate in the determination 
of damping, and in error by an order of magnitude at two octaves away 


from the resonant frequency. 


Fig. 3.16 shows the rough computed spectra of the pitch signal whose 
auto-correlation function was given in Fig. 3.7, А good agreement with the 
exact transform of the approximate damped cosine wave is demonstrated, 


except for the erratic behavior at low frequencies. 


The above examples were fairly straight-forward because only a 
magnitude function had to be approximated. In the case of a cross-power 
density spectra, the difficulty is increased by the need to match a phase 
curve as well as one of magnitude. This is identical to the problem in con- 
trol systems design, when a transfer function is to be found for a system 
described by experimentally determined magnitude and phase curves of 


Sinusoidal response. 


This procedure is strictly one of intelligent trial and error. Several 
good references are Chestnut and M crum who treat the basic problem, 


20 
and Truxal , who discusses more refined estimates in the iterative cycle. 


A plot of magnitude and phase for a single pole is given in Fig. 3.17, 
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and for complex poles in Figs. 3.18 and 3.19. 


In the case of cross power spectra, one desires to determine the 
LHP and RHP poles and zeros. The following table lists the correct sign 
and magnitude of the logarithmic magnitudes and phase angles of Figs.3.17-19. 


LOG MAGNITUDE PHASE 
RHP ZERO Е > 
LHP POLE + + 
RHP ZERO = + 
RHP POLE t - 


Table 3.1. Sign of magnitude and phase terms 


If the cross power spectrum Ф. =) is to be found, it will usually 
be true that its positive or LHP poles will correspond to those of Ф ll 
and its negative or RHP poles will be those of б. AS). 


The roll-pitch cross-correlation functions shown in Figs. 3.8 and 
3.9 have a digitally-computed cross spectra given in Fig. 3.20. An appro- 
ximation over the range of validity appears superimposed. To achieve the 
proper phase response, it was necessary to introduce a RHP zero and LHP 
pole pair which were complex and mirror symmetric. Consideration of 
Table 3.1 indicates that the magnitude effect of these terms is zero, but 
that the phase effect is additive, giving a net phase shift of -360 degrees 


at high frequencies. 


The analysis of cross-spectra with logarithmic plots of frequency 
behavior can make no use of the co- and quadrature spectrum records. 
Accordingly, it is urged that magnitude and phase plots become the standard 


presentation. 


3.4 Conclusions 


(1) A review has been made of the normal methods used for the approxima- 
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tion of correlation functions and power spectra density curves with terms 
usable in a linear theory. No distinct preference has been recorded for 
either method -- in fact, they supplement each other and in the determina- 
tion of the dominant frequency and damping values both are equally efficient. 
(2) The method currently advocated for the computation of spectra 

from sampled correlation functions causes considerable error which could 
be lessened with a more sophisticated computation scheme. 

(3) The use of a pole-zero approximation to observed wave-height spectra 
appears much more desirable, especially if the damping ratio Y remains 
sensibly constant in the dominant frequency C than is the preoccupa- 
tion with the empirically-derived Neumann spectrum. 

(4) Statistical analysis of ship motion does not generally reveal useful 
knowledge of the input wave forces and moments because of the high degree 


oí filtering which the ship system performs on this signal. 


Inge 





CHAPTER IV. 


CONTROL OF SHIP MOTION 


4.1 Introduction to multi-dimensional ship control 


Let us consider a hypothetical ship which is equipped with individually 
movable fins at both the bow and stern and a rudder. Motion sensors are 
available to measure any or all of the degrees of freedom of the ship. The 
design sea is rough, and the peak wave forces exceed the maximum forces 
from the hydrodynamic appendages because of the latters' finite capacity 
and breakdown at large angles. It is desired to minimize the average 
weighted sum of total squared displacement error and total squared velocity 
of a particular point aboard ship in three dimensional motion, perhaps at 
a critical missile launching station, by simultaneously controlling each of 
the five control surfaces. The interactions in the various degrees of free- 
dom of the ship are to be accounted for, as well as the statistical interrela- 
tionships of the seaway. In order to maximize the stabilization per dollar 
expended for the appendage area, these fins and rudders are to be driven 


full throw. 


In short, this problem is one of the most difficult in control engineer- 
ing. This chapter will present the core of a theory which will allow solution 
of this and other related stabilization problems in a near-optimum and 
easily-mechanized fashion. The theory will be applicable to surface ships, 
submarines, hydrofoil craft, "hydro-skimmers" -- in short, in any marine 
vehicle problem where the sea exerts the disturbing random force and it 


is desired to stabilize with several control variables. 


4.2 Inadequacy of a linear theory 


Chapter 2 has dealt with random processes and the optimum systems 
to operate on them. It is logical to inquire if there is an optimum linear 


System for the ship control problem. In a related example, Fig. 4.1 shows 
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an input random process v, a Scalar fixed system H(s), and a controller 
KC(s) which acts to reduce the output error e. In this simple case it is 


obvious that including an infinite gain in the controller, (ie: K > °S ) 





Fig. 4.1 System subject to disturbance 


will cause the output or error to be zero (if the resulting system is stable) 
since 


E(s) _ H(s) 1 as K — e 
V(s) ` 1 +K C(s) H(s) 77? K C (s) 





Unfortunately, in the ship motion problem infinite gain is impossible 
in the "controller", whose output is viewed as the force of the fins or rudder. 
This infinite gain would imply an infinite force output for a finite error, which 
is patently not the case. The real limitation for ship motion control is the 


peak force which can be exerted by the hydro-dynamic control surfaces. 


From Nen a theory exists which gives the optimum linear 
controller when its output is constrained to have a certain average power 
output. This appears superficially to be a valid approach, but actually it 
is forcing the physical problem to fit the theory, rather than vice-versa. 

It further limits the corrective force available instead of making intelligent 


use of the known limitation. 


The designer of a ship pays a considerable penalty in space and 
weight for large fins and rudders and their associated high-power actuating 
devices. Because of this high price, he naturally desires to maximize the 
stabilization obtained from his installation, even at the expense of using 
quite sophisticated signal processing equipment, including perhaps many 


motion sensors and a digital computer. With this assumption, we must 
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discard our hopes for a completely linear theory since experience has 
ses that non-linear control systems in general have much better 
performance in practice than linear systems, if the former are adequate- 


ly designed. 


If an error or the difference between the actual state of the ship 
and the desired equilibrium position exists, the control system will attempt 
to reduce this error by exerting corrective effort. It is intuitively appeal- 


' at every instant if its 


ing that this corrective action would be ''better' 
magnitude were greater. That is, a maximum-effect, full-throw or "bang- 
bang" control system, which at every instant uses the maximum lift force 
possible at every hydrodynamic surface, will be postulated in the develop- 
ment to follow. A linear mode of operation in an arbitrarily small region 
of system error will be permissible, however. This full-throw type of 


control has been shown to be optimum in a class of problem which is con- 


cerned with the settling of a free system to equilibrium in minimum time. 


4.3 A multi-dimensional approach for saturating controls 


The design philosophy to be presented in this section was first 
given in Ref. 1, and is basically an extension of the so-called "Second 
Method of Lyapunov", avery general theory of dynamic behavior of linear 
and non-linear systems which is currently receiving considerable atten- 
tion in the ие E. A summary of this approach is presented first, 


followed by a more complete development of the theory. 


In general, we first attack this problem by specifying a quadratic 
criterion of instantaneous error for the ship, which could be, for example, 
the weighted sum of all the various squared displacements and velocities 
(ie: state variables) existing at a certain time. More formally we postu- 
late an instantaneous measure of error with the quadratic forie zn Qu 


where x is the state vector of the ship from Section 2.2, and Q is a non- 


negative definite symmetric matrix. 
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In control, we are interested in the future as well as the present, 
and a measure of the future error of our instantaneous state is quite 
naturally given by the integral of the instantaneous error criterion over 
all future time. Since a random process is exciting the system, the ex- 
pected value of future components of x is used in place of a determinate 
x. Also, future values of x are a function of the random process state 
variables, as well as those of the ship. By considering y to be the state 
vector of a linear system which includes both the ship and the random 
process generating model, the Second Method of Lyapunov is used to de- 


termine a symmetric matrix P, such that 
oO 


T T 
у Ру = Pas Q x 
=, 


Thus, at every instant oftime, a single number ae is a measure 
of all knowable future error, neglecting future values of control inputs. 
The various control signals are based on the requirement that the negative 
rate of change of this measure of future error be maximized continuously. 


The advantages of this approach are: 


(1) The theory becomes more valid as the sea becomes rougher 
and the saturated outputs of the force transducers become relatively feeble 
with respect to the wave excitation. 

(2) The design procedure is simple to execute and resultsina 
completely linear system except for the saturation constraint. 

(3) The resulting system is guaranteed to be stable, 

(4) The multi-dimensional and inter-coupled aspects of the problem 


do not appreciably complicate the solution. 
Now, to consider in detail the aspects of this procedure. 


First, it is necessary to select the instantaneous criterion of error, 
O This is the area in which the skill and knowledge of the designer 
is applied with best effect. As in all optimization procedures, the critical 
and often least understood part of the design is deciding what is really 


wanted. 
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One rational choice is for Q to be a diagonal matrix with positive 
elements, 2r ; ау would then be 2 А. Б. The Zw. could be selected 
2a 


so that x is approximately proportional to the total energy of the ship. 


Since kinetic energy is represented by = M y? and potential energy by 


2 
1 
2 K x^, with K an effective spring constant, the ^ weighting terms could 


be selected proportional to the acceleration and displacement coefficients 
in the differential equations relating the velocity and displacement state 


variables, respectively, neglecting the cross-coupling. 


Another potentially useful selection of the matrix Q would represent 
N 
E Vo: M ; 
the weighted value of squared velocity inthree dimensions, at a single 


error by a x, + the squared disturbance displacement plus 
point "a" on the ship. To obtain this measure easily, suppose that the 
vector F locates the point a from the center of gravity of the ship. If the 
vector v represents the instantaneous disturbance velocity of the center 
of gravity -- that is, the actual velocity less the mean forward velocity, 
the vector Ор the angular displacement (for example, the roll 
displacement ф would be a vector along x axis with a magnitude equal їо 
ф in radians), the vector x is the disturbance displacement of the center 
of gravity, and the vector Wis the total angular velocity of the ship, then 


the motion of the point a is given by 


en — > o> 
х = x t Oxr 
a 

SC Z5 — „мм 
ЕЕ. Ox Гг 


for small angles of displacement. 


Let the state of the ship be defined as in Section 2.2, 


x, = Roll velocity, p Xo 7 Surge velocity, u 

X) = Yaw velocity, r Xp = Roll displacement, 

xq = Sway velocity, v ха = Yaw displacement, ل‎ 
x, = Heave velocity, w Хо = Heave displacement, 2 
х = Pitch velocity, q Xio” Pitch displacement, © 


Then, from evaluation of the vector equations, 


os 
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where ro is p locate the point of the ship to be stabilized from the ship 


center of gravity. 


If 
X - R x | 
a = 
Then, 
ZX 2 = m = ес 
A 


and 





x ha -h Yu -ҮхҒу O | 
- hy + f © -Yafa =“ | 
rx | O O о | С) 
O O | EJ. cS 
“Yuta О -ҒҰ К+» re | 
BEL ye Pe Or: COC BE EH DE о 
| YE -Ae y -RA 
O Bea +R o — (> 
| ry O | - fx 
| -RA -UR -R WR 


Having determined, by any means, a measure of instantaneous 
error, x Qx, it is now necessary to evaluate the integral of this quantity 
over all future time. To do this, we use particular results from the Second 
Method of moe First of all, let us define the system under con- 
sideration, Fig. 4.2 illustrates the creation of the random forces and 
moments of the sea from a white noise driven model G(s). These forces 


E | 
excite the ship, H (s), anda controller, sensing the motion, adjusts its 
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various control variables with full throw. 
Section 2.2 showed how the matrix differential equation 


E х = Ах жр. Ес + р. Е 

was formed for the ship alone. If the concept of system is enlarged to 
include G(s), which of course has its own state variables, then the differ- 
ential equation of this new system can be written as 


d 
=> = 1 

ar у By DE (4. 1) 
where y now contains, besides the ship state variables, the state variables 
of G(s). F is the excitation vector, having controlled force inputs and 


white noise as its components. 
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Fig. 4.2 The multi-dimensional control problem 


If y has for its first ten variables the corresponding state variables 
of x, which is certainly allowable, the instantaneous error criterion ON 
can be formed by using the Q matrix discussed earlier in the first ten rows 
and columns, and filling in the remainder of the matrix with zeros. This 
is equivalent to saying that we assign no instantaneous error to a state 
variable of the sea, since we cannot control it. 


If the relationship 


T T 
f » Qydt- y Ру 


A 


= 





is to be satisfied, then by differentiating the above equation we obtain 
T d [f т ] d т) ma ) 
= = ашаны = ane + P nn 
05 УУ EY c ЖАР) dt У 
ex 


T | 
Now, in the computation of J y Qydt we are concerned with 


the expected value of future y, disregarding future control action. Hence 


da Y^ BY 
Substituting, 
Е во ep, 
Therefore, 
Bie pte Pee = - © 


The elements of the symmetrical matrix P can be solved for 
n(n + 1) 
2 
equations, where nis the dimension of y. 


numerically by the independent equatibns in the above matrix 


The criterion to be used to fix the sign of the control variable is 
that the measure of future error, y н P y, is to have the maximum possible 


rate of decrease at every instant. In other words, 


ы 
KENT mo Ede т ту чу pld} 
dt | Pyf- E y j M PF = а) РЧ 


is to have a maximum negative rate of change. But 


ау . 
dt - E SENDER 
Substituting, 


d T Terr Por 
E ру} 2 {утв ter рт} ру 


The effect of the control variables should be such that ЕТ рт Р y! 
is as negative as possible. This will obviously occur when F., the Ка con- 
trol variable, is at its maximum, and has the opposite sign of the ү) com- 


ponent of р! Py. Thatis, 


I ae IF, max | sgn {oT py J, (4. 2) 


Si 





where sgn x = +1, Xx <O 


= -], x <O 
Eq. 4.2 is the desired control law. 


T ee 
The numerical coefficients of - D P provides a weighting of the 
measured values of displacement and velocity in the various degrees of 
freedom of the ship, and also of the state variables of the random process 


generating model, G(s). 


The measurement of these latter state variables poses an interesting 
problem. First of all, with the use of accelerometers, rate gyros, etc., 
it is feasible to measure the motion of the ship. Thus, the equations of 
motion of the ship, summarized by H(s), allow direct reproduction of 
the wave forces, F y by suitable weighting of the outputs of each of the 


motion sensors. 


This six-dimensional random process will have been previously 
analyzed by statistical methods to obtain a matrix of power density spectra. 
Suppose that there is a small but finite noise in each of the signals re- 
presenting the wave forces, which would certainly be the case in practice, 
and an approximation to the statistical nature of this noise is available 


(probably a white noise assumption) . 


Even if there are steady state sinusoidal components of noise, 
which might occur with sustained vibration of the ship's structure, for 
example, this approach is still applicable and results in filters which 
tune out the offending sinusoidal components. By representing the input 


spectra of the wave forces as measured by Ф (s) = ф (s) J Ф (s), 








where @..(s) represents the signal or force spectra, and @ (s) the un- 


an 


wanted noise spectra, then Q As) can be factored into G(-s) G (s) by 


the methods of Ref. 1. 











G(s) is separated, by partial fraction expansion, into S(s) and N(s) 


terms respectively, where for the assumption of uncorrelated white noise, 
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N(s) willbe a numerical matrix. Then, the multi-dimensional feedback 


system shown in Figure 4.3 is constructed. 


STATE VARIABLES 


Y to ۷ 
W COS E f 1 з: 
Noisy - . OF 

WAVE N ($) > (5) ° WAVE FORCES 
_ РОВСЕ $ 

У, + ый 
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Fig. 4.3. Optimum multi-dimensional filter used to produce state 
variables. 
This shows the configuration of an optimum multi-dimensional 
filter, as was discussed in Sections 2.3 and 2. 4. In this application, the 
output or filtered wave force signals are not used, but the inner signals 


of S(s) are the desired state variables. Fig. 4.4 outlines the proposed 


scheme of control. 
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Fig. 4.4. A full-throw control system for multi-dimensional ship 
stabilization 
Here, K, weights the outputs of the motion sensors according to 
the equations of motion of the ship to calculate the external forces. N 
and S act as an optimum filter to reproduce the state variables of the wave 


force random process. Ко weights the state variables according to Eq. 4.2, 
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with the matrix, - DP: 


ERR has conducted an analog computer investigation of a 
particular scalar problem, stabilization of a second order system, using 
this theory. His results show a superiority of this approach as the system 
damping ratio becomes low and the ratio of maximum control force to 
r.m.S. disturbance force becomes small, compared with other "optimum" 
design procedures. As this domain of superiority is especially applicable 
to ship motion, considerable promise is offered for the multi-dimensional 


Case. 


Thus, to summarize, the multi-dimensional full-throw controller 
is made up of signals which it is practical to measure and performs 
essentially a linear computation, except for the binary nature of the output 
control commands. A fair amount of analog measurement and calculation 
is called for, but the basic problem itself is a very complex one, perhaps 


among the most difficult encountered in control systems engineering. 


4.4 Some problems of practical multi-dimensional ship motion control 


The design presented in the previous section is in a sense a central- 
ized one -- it considers the ship as a whole and makes continuous decisions 
as to how each control variable should best be switched. In this section 
we will consider several problems which arise in the practical mechaniza- 


tion of this approach: 


(1) Practical force transducers, such as hydraulically driven fins 
and rudders, require a finite time to reverse the direction of maximum 


force and cannot "switch" instantaneously. 


(2) The actual hydrodynamic force exerted on a fin depends con- 
siderably on the relative direction of motion of the local water flow and 


not entirely on its mechanical displacement with respect to the ship. 


(3) Practical motion sensors can be expected to deliver a signal 
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which is only an approximation to the true ship motion because of internal 


noise, cross-coupling effects, and local vibration levels. 


In considering the first objection, we assume that the actuating 
equipment is sufficiently powered to allow full travel in a time equal to 
or less than the mean switching interval. If this assumption is not satis- 
fied, then the system is rate limited rather than amplitude limited and 
the design premises are shaky indeed. In this case, certainly the capacity 
(ie: area) of the fins would not be efficiently used and either a reduction 
in their size or an increase in the power rating of the fin mover would 


appear to offer a more economically optimum design choice. 


The output of the weighting operation (n! P yf. it will be recalled 
from Section 4,3, is the multiplier of the e control force which determines 
the rate of change of the error criterion. When it has the value zero, ob- 
viously the force exerted at that instant contributes nothing to the reduc- 
tion oí this criterion and might as well be zero. That is, since a finite 
time is required to go from full up to full down, a near optimum choice 
of when to commence the downward travel would be so as to make the 
force pass through zero at the same time as [- DE P y}; does. This would 
minimize the effect of the necessary small values of force required in 
a finite-time passage through zero by weighting them with small values 


of multiplier. 


Thus, the signal for the ie force variable to commence changing 
sign should be the prediction of sgn [- De P У} at T: seconds in the future, 
where 2 T: is rc for full XN of thei fin. This will be 
given by sgn [. D Pe 1 у}, where e i is the matrix exponential 
defined in Section 2.2, and Bis given in Eq. 4.1. The switching computer 
K, shown in Fig. 4.4 would then have an sue row of - рїр е S : 

The second problem, that of relative flow effects, is not so easily 
handled. The trouble is pictured in Fig. 4.5, which shows the lift break- 


down of a typical hydrodynamic surface. Suppose that a feedback signal 
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is available of the actual lift on the fin, which is customary in contemporary 


fin installations. If the maximum force available and commanded ıs B, and 





LIFT FORCE 





DISPLACEMENT RELATIVE To 
INSTANTANEOUS FLOW 


Fig. 4.5. Ambiguous force-displacement relation in fin breakdown 


the actual force measured is given by the dashed line in Fig. 4.5, then 
it is not clear whether the relative displacement is at Aor C, and the 
direction for corrective motion is unknown without (for example) additional 


pressure distribution measurements, 


This problem is quite difficult, and typical of a general "peak-hold- 
ing'' class of control problem. While the mechanization of the solution to 
this interesting dilemma will not be elaborated upon further, one possibility 
would employ a small high-frequency sinusoidal ''dither' on the fin -- or 
on a lift flap -- and the phase of the resulting sinusoidal measured force 
would distinguish between the case of point A and C. If the measured al- 
ternating force is in phase with the dither, then fin breakdown has not 
occurred, and, in fact, this sinusoidal force signal represents a valid 
error signal for a positioning control system for the fin since it goes to 


zero at the desired peak output. 


The compensation of the noisy outputs of motion sensors is a 
standard problem of instrumentation engineering. If the noise can be 
statistically approximated by a power spectra, the theory in the previous 
section for optimum recovery of signal information is immediately appli- 
cable. Also, a convenient help in this respect would be the installation 
of sharply tuned filters at the instrument outputs which would reject sinus- 


oidal components at dominant hull vibration frequencies excited by slamm- 


zo 





ing, for example. 


4,5 Adaptive measurement of wave force Spectra 


The key assumption in the previous theory is that the random process 
which represents the six-dimensional wave forces on the ship is stationary. 
This is, of course, not true, since ship course and speed changes vary 
the spectra as felt by the ship and since the properties of the sea change 
over a period of hours. However, the time variation of the parameters 
of the sea is very slow, so that the process can be termed a quasi-stationary 
random process and the previous theory is valid, providing that the current 


statistics of the sea are known. 
There are three possible approaches to this problem: 


(1) Since the statistics of the sea are not known, neglect them and 
their associated control theory, and carry out designs by more conventional 


(if less effective) means. 


(2) Consider the worst case, statistically, which the ship could 
encounter (that is, excitation at the resonances of the ship) and design 
for this condition, accepting the degradation of performance at less severe 


sea states, 


(3) Determine the dominant characteristics of the sea spectra, 
and attempt to measure them continuously and alter the control system 


to reflect these changes. 


Some possibilities of the third approach are to be discussed in 


this section. 


It is difficult to discuss intelligently the properties of the wave- 
force random process without experimental verification. However, there 
have been some statistical analyses of the sea which indicate considerable 


promise in adaptive measurement of wave forces, 


. 12 À : 
Voznessensky and Firsoff made an extensive series of measure- 


E ee 





ments of wave heights at sea using a wave pole, and performed auto-correla- 
tion analyses. They found that the correlation function was strikingly ap- 
proximated in all cases by a damped cosine wave, к, enu: DI . From 

the discussions of Chapter II, it is clear that a generating model with 


the transfer function 


K, 
s^ Ta ST $ 
n n 
excited by white noise -- a standard second order system -- could be used 


to reproduce the process. The important result of these measurements 
was that the ratio of = remained essentially constant at .21, despite the 
wave-heights and directionality of the particular wave system measured. 
The constant ratio mer ‚ 15 the non-dimensional damping ratio 5 
in the standard E TÉS constancy appears to be a natural property of 


waves at sea, and tremendously simplifies the measurement of the spectra, 


as will be discussed. 


The wave spectrum, if V is constant, is totally specified if the mean- 
square value of the signal and a is known. The mean-square value of 
the signal is easily determined by a squaring and filtering operation of 


an electrical signal representing the process. 


The M value can be found through a process of zero-counting, or 
determining the average number of zero-crossings of the signal per unit 
time. From the work of Be Ehrenfeld, et a and ner it 
can be shown that for a standard second-order system excited by white 


noise, the expected number of zero-crossings per secondis “п 
7 
The possibilities of evaluating the dominant characteristics of a 


random process by squaring and averaging, and zero-counting, is a power- 
ful tool for continuously evaluating the parameters of the sea. It is hypoth- 
esized that the dominant quantities of the multi-dimensional wave force 


spectra, since they arise from wave motion, can be measured by these 
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simple operations. As a simplified example, the effective frequency of 
excitation of the pitch-heave system, and that of the roll excitation, could 
be measured continuously by zero-counting the calculated forces as re- 


produced in Fig. 4.4. 


The ability to make good statistical estimates of the parameters 

of the input random process without the need for a full-scale correlation 

or spectral analysis is a considerable boon to the control systems designer. 
It is felt that a valid approximation to the third alternative at the first of 
this section would be, on the basis of future experimental measurements 

of spectra at sea, to fix many of the parameters as sensibly constant and 

to measure some of the most important characteristics -- in particular, 
"frequencies" of excitation, мше and design a control system which would 
be nearly optimum in the sense of this chapter and which would be contin- 


uously "adapted" to changes in ship's course and speed and sea spectra. 


4.6 Summary 


The ship motion control problem is one of great complexity because 
of the inter-coupling of ship motion, the correlation between the various 
wave forces, and the saturation of control transducers. This chapter has 
presented one solution to this problem which has the virtues of near optimal 
performance, and the possibility of mechanization with a reasonable amount 
of computing equipment. Needless to say, only the germ of the idea appears 
here, and the real test of its validity must await accumulated measurements 
of ship force spectra and subsequent analog computer simulation of the con- 


trol system. 
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СНАРТЕВ У. 


CONCLUSIONS AND RECOMMENDATIONS 


5.1 Introduction 


The penalty for making a frontal attack on such a complex problem 
as ship motion in a random seaway is necessarily a compromise between 
a desire to examine all facets in detail and a need to communicate the es- 
sential ideas clearly. This chapter will first attempt to outline the major 
viewpoints which have comprised the backbone of this thesis. Then, a 
specific program of proposed research to implement these ideas will be 


presented. 


5.2 Major ideas 


The basic idea which threads through all past chapters is that the 
random wave forces and moments on a ship at sea are the dominant quan- 
tities of interest in any study of the complicated behavior of multi-dimens- 
ional ship motion. With the development of suitable theory in Ref. 1 to 
analyze the statistical interrelations among these forces and moments, 
all useful information regarding (1) average motion levels , (2) short term 
motion prediction, and (3) optimum control decisions can be found. ш short, 
a path is now opened which will explain all that can be explained about ship 
motion within the limits of current mathematical knowledge and, what is 
perhaps more significant, the needed measurements can be easily made 


with standard instrumentation and data reduction methods. 


Chapter Il has briefly shown that the ship can be described in the 
two basic models for linear systems -- the matrix transfer function in 
Laplace transform notation and the matrix differential equation, the latter 
using the significant concept of "state". The idea of a white noise driven 


model has been introduced, and the statistical measures of power density 
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spectra have been interpreted as means for determining this linear model. 


The practical problem of approximating experimental spectra with 
rational Laplace transforms has been examined in Chapter III. Besides 
presenting some graphical aids in this task, it has been demonstrated 
that methods now in use for statistical data reduction using a digital com- 


puter introduce unnecessary error in the final spectra. 


Also, it was shown that a considerable degree of statistical correla- 
tion exists between the roll and pitch motion of a ship at sea, the major 
portion of which is hypothesized to arise from the correlation between the 


respective wave excitation forces. 


Chapter IV has developed in some detail a proposed scheme of 
simultaneous control of many force transducers to stabilize a ship at sea. 
With the virtues of near-optimality, full-throw control, case of instrumen- 
tation, stability, and consideration of all interactions in the ship motion 
and the input forces, this scheme offers considerable promise of extract- 


ing near the maximum in stabilization from fins and rudders. 


Preliminary analog computer investigation reported ы 
gives substantial support to this design concept for ship control. The 
peculiar problems of relative water velocity effects, hydraulic rate satura- 


tion, and adaptive spectral measurement have been discussed. 


9.3 A proposed research program 


The following is a proposed program of needed and unified basic 


research in the areas of ship behavior at sea and ship motion control. 


A. Ship parameter measurement 


Determine, experimentally and analytically, the significant co- 
efficients of the linear differential equations (as a function of forward 


velocity) which govern the total motion of the particular ship type under 


So 








investigation. Use free or constrained models with a suggested greater 
use of initial condition response. Also, statistical analysis of measured 


ship motion at sea will yield useful verification of the small-scale results. 


B. Wave force measurement 


On a ship at sea, or on a free model in a random wave generating 
facility, make sufficient weighted motion measurements to infer -- through 
the assumed linear differential equations of motion -- the wave force or 
moment excitation in each of the six degrees of freedom. An alternate 
experiment could make use of a completely rigid model in a random wave 


tank, with the desired forces measured directly. 


C. Immediate wave force analysis 


From the six-fold record of wave forces, conduct a complete cross- 
spectral analysis, either by digital computation of correlation functions 
and subsequent transformation, or by direct analog measurement. Approxi- 
mate these spectra by rational transforms. Factor, preferably with a 
digital computer, the matrix of power density spectra with the method 
introduced in Ref. 1. The result, a 6x6 matrix of transfer functions, 
completely characterizes the random nature of the sea and can be imme- 


diately used in analog computer studies. 


D. Extended wave force analysis 


Collect a number of matrix descriptions of the sea. Attempt to 
determine relatively invariant parameters, as well as the most significant 
varying ones. Define a typical or design sea, possibly on a worst-case 


basis for the ship under consideration. 


E. Relation to wave-height directional spectra 


Conduct theoretical and experimental analyses of the relation 
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between the observed statistical properties of wave forces on a ship or 


model and the directional spectra which characterizes the sea. 


F. Multi-dimensional control 


Using the matrix model developed under steps C and D above, con- 


duct analog computer investigations with the following objectives: 


(1) To evaluate the control theory outlined in this thesis, as well 
as to stimulate other approaches to the multi-dimensional stabilization 


problem. 


(2) To determine the importance of using knowledge of the input 
random process as an aid to control, as well as the penalties associated 


with imperfect or no knowledge of it. 


(3) To determine the optimum balance between large fast power 
actuators of rudders and fins and large surface areas of these appendages, 
and also the reduction in stabilization because of rate saturation and lift 


breakdown. 


G. Ship motion analysis 


From accumulated knowledge of wave-force spectra, determine 
the significant parameters of the ship which affect seaworthiness and 


emphasize undesirable motion. 


H. Adaptive measurement of seas 


Conduct analytical and experimental investigations into adaptive 
measurements or "tracking" of significant parameters of random processes, 
with the goal of providing a shipboard control system with a self-adjusting 


capability as sea conditions change. 


5.4 Conclusions 


The above program, which unifies the central themes of this thesis, 
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is believed to offer a promising and practical roadmap into the vital areas 
of ship motion and control. Many problems and an unbelievable amount 
of data remain to be analyzed and digested, but no valid approach to such 


a complex problem as this can avoid such a large-scope investigation. 
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